Energy spectra stemming from interactions of Alfven waves and turbulent eddies 
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We present a numerical analysis of an incompressible decaying magnetohydrodynamic turbulence 
run on a grid of 1536 3 points. The Taylor Reynolds number at the maximum of dissipation is 
~ 1100, and the initial condition is a superposition of large scale ABC flows and random noise at 
small scales, with no uniform magnetic field. The initial kinetic and magnetic energies are equal, 
with negligible correlation. The resulting energy spectrum is a combination of two components, each 
moderately resolved. Isotropy obtains in the large scales, with a spectral law compatible with the 
Iroshnikov-Kraichnan theory stemming from the weakening of nonlinear interactions due to Alfven 
waves; scaling of structure functions confirms the non-Kolmogorovian nature of the flow in this 



range. At small scales, weak turbulence emerges with a k A 
referring to the local quasi-uniform magnetic field. 



spectrum, the perpendicular direction 
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Magnetic fields permeate the universe, and with the 
increased resolving capabilities of instruments, a flurry 
of details on highly complex flows emerge. The origin 
of such magnetic fields, the dynamo problem, is still 
not fully understood, in particular when the magnetic 
Prandtl number Pm - the ratio of kinematic viscosity 
v to magnetic diffusivity r\ - differs substantially from 
unity, as it does in the interstellar medium (Pm 3> 1) or 
in convective regions of stars and liquid cores of planets 
(Pm -C 1). Numerous theoretical and numerical studies 
have been written (see [ljj), and recent laboratory exper- 
iments address this fundamental problem as well. Here, 
we take a different approach by assuming that a dynamo 
mechanism works, presumably leading to approximate 
cquipartition of the magnetic and kinetic energies, as ob- 
served for example in the solar wind. We then ask ques- 
tions about the nature of the nonlinear dynamics of such 
a flow. The complete problem (a dynamo taken all the 
way to the nonlinear stage at a high Reynolds number 
as encountered in nature, and for Pm differing substan- 
tially from unity) is still out of reach, but it might be 
among the key outcomes of petascale computing facili- 
ties and experimental devices currently being developed. 
We specifically address in this paper the nature of the 
energy spectra that arise in such conditions. A more de- 
tailed account of the development and evolution (up to 
the peak of enstrophy) of structures and correlations of 
the flow will be presented elsewhere. 

The magnetohydrodynamic (MHD) equations read: 



<9v „ 1 „ m 
— + v Vv = VV 

at po 



j x b 



^V 2 v, 



9b 

— = V x (v x b) +7?V 2 b; 
at 

v is the velocity, b is the magnetic field, j = 
the current density, V is the pressure, po = 
sity, v = 7} = 2 x 10~ 4 , and V • v = V • b 
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= V x b is 
1 the den- 
= 0. One 



can also write these equations in terms of the Elsasser 
variables = v ± b, of energy E^ and flux e ± . We 
solve Eqs. (Q] [2]) in a three-dimensional box of length 
Lq = 2?r using periodic boundary conditions and a pseu- 
dospectral method, dealiascd by the standard 2/3 rule; 
minimum and maximum wavenumbers are k m i n = 1 and 
kmax = AT 1/3 /3, with N = 1536 3 grid points. At all times 
ko/kmax < 1) where ko is the dissipation wavenumbcr. 
The initial conditions for the velocity and magnetic fields 
are constructed from a superposition of three Beltrami 
(helical) flows to which smaller-scale random fluctuations 
are added with initial kinetic and magnetic energy Ey = 
Em = 0.5 (with respective spectra Ey(k) and £j\f(fc)); 
magnetic helicity Hm = (a • b) « 0.45 (b = V x a where 
a is the vector potential, and the brackets denote vol- 
ume average), and cos(v, b) = (v • b) ( | v 1 1 1» | ) 1 as 10~ 4 
(see Q for details). The computation is stopped when 
the growth of the total dissipation saturates (t = 3.7), at 
which time the Reynolds number based on the mechani- 
cal integral scale is R e = ULy jv as 9200, and that based 
on the mechanical Taylor scale is R\ = UAy/v as 1100; 
U is the r.m.s. velocity, the integral scales are de- 
fined as Li = 2nE~ 1 j k^ 1 Ei(k)dk, the Taylor scales 
are A; = 2n(Ei/ J tfE.^dk) 1 ' 2 , and i is either V or 
M. These scales at t = 3.7 are Ly sa 2.6, Ay as 0.31, 
Lm ~ 3.1, and Am ~ 0.39. It was shown in [2| that at 
early times, current and vorticity sheets form, and fur- 
ther roll-up, fold, or pile-up. Here we focus on the fully 
developed regime close to the peak of enstrophy. 

The total energy flux H(k) as a function of wavenum- 
ber at t = 3.7 is shown in Fig. Q] (inset), together with 
the total, kinetic, and magnetic energy spectra (E(k) = 
Ey(k) + E M (k)), compensated by either fc^ 5 / 3 (Kol- 
mogorov spectrum,hereafter K41) or fc~ 3 / 2 (Iroshnikov- 
Kraichnan, or IK [3j); the latter takes into account the 
slowing down of nonlinear transfer because of Alfven 
waves. The scaling of E(k) and Em (k) is close to fc -3 / 2 in 
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FIG. 1: Energy spectra E(k) (solid), Em (k) (dash) and 
Ey(k) (dot) compensated by fc _3//2 ; E(k) compensated by 
fc~ 5//3 is also shown (dash-dot). In the inset is the total en- 
ergy flux Il(fc), indicating the extent of the inertial range. 



a range of scales covering almost a decade. Energy spec- 
tra in MHD have been reported in the literature, both 
for decaying and forced flows, for data stemming from 
either numerical simulations or observations (mostly in 
the solar wind). In most cases, a spectrum close to K41 
is found, although there are times when the IK solution is 
preferred Q , and in some cases both scalings may be ob- 
served Q. Moreover, it is not clear whether universality 
obtains (in which case, for very large Reynolds numbers it 
would either be K41 or IK), or whether there are classes 
of universality, the boundaries to which are not neces- 
sarily known (see, e.g., the discussion in Q). Note that 
different spectra have already been observed for differ- 
ent forcings in reduced MHD (RMHD, an approximation 
valid in the presence of a strong magnetic field) Q, and 
recently also reported for MHD runs with a strong im- 
posed magnetic field. Note that the flow in our run has 
the largest R\ obtained in a direct numerical simulation, 
at the expense, though, of not being stationary. 

To clarify the issue of which spectra should arise in 
MHD turbulence, at least two things can be done. On 
one hand, one computes isotropic high-order structure 
functions, in which case the differentiation between K41 
and IK is enhanced, as already found in Q. On the other 
hand, one can quantify the degree of anisotropy in the 
flow, something that will be discussed later. We define 
the longitudinal (parallel to the displacement 1) struc- 
ture function of order p for the quantity f as Sp(£) = 
([Sf L (l)] p ) = ({[f (x) - f(x + 1)] • 1/£} P ). If the field is 
self-similar we expect a scaling Sp(i) ~ in the iner- 
tial range, where ^ are the scaling exponents. In hy- 
drodynamic turbulence, K41 theory predicts ( p = p/3, 
and in MHD the IK theory leads to C, p = p/A. In prac- 
tice, these relations are modified by intermittency cor- 
rections due to extreme events at the small scales. The 



FIG. 2: Scaling exponents C* for z ± = v±b. The K41 (dot) 
and IK (dash) predictions are also shown. Note ~ 1. 



exact scaling laws for MHD turbulence derived in [9j], 
(fe^l^z^l)! 2 ) = -4e ± ^/3, can be used to define the 
inertial range and to improve the estimate of the scal- 
ing exponents, as is often done in hydrodynamics when 
the extended self-similarity (ESS) hypothesis is used [Toj ] . 
These scaling laws, and structure functions up to order 8, 
were computed at t = 3.7 for increments in the x, y, and 
z directions, and averaged. Figure O shows the resulting 
scaling exponents for the Elsasser fields z ± , computed 
in the range in which the exact scaling laws hold. 

Consistent with the computation of the energy spectra, 
the scaling of z ± is closer to the IK prediction, with = 
0.831 ± 0.006, C 3 " = 0.806 ± 0.005, Ct = 0.985 ± 0.008, 



C7 = 0.958 ± 0.007. As stated in [9j, the exact scaling 
laws in MHD involve correlations between the velocity 
and magnetic fields (or equivalently between z^), and 
such correlations play a role in the dynamics of the flow. 
In particular, they can induce a shift from a standard K41 
coupling; this effect has been modeled in ll| as a scale 
variation of the angle between the velocity and magnetic 
field, i.e. in the relative cross helicity cos(v, b). From 
the values of and , the possibility of a Kolmogorov 
scaling hidden by a bottleneck (as reported in 12|) can 
be ruled out for this simulation. However, even when the 
scaling exponents in the inertial range for p = 3 and 4 are 
closer to the IK prediction (with strong corrections due to 
intermittency), it is worth mentioning that the structure 
functions after using ESS do not show the same scaling at 
all scales. This is different from the hydrodynamic case 
where, after using the ESS hypothesis, a unique scaling 
at all scales is observed, even in the dissipative range. 
This brings us to our second point, linked to anisotropy. 

Indeed, it is known that in MHD anisotropy develops 
at small scales under the influence of a large-scale mag- 
netic field. Several models have been written to explain 
how this might affect the small-scale dynamics. The weak 
turbulence (WT) theory [l3| leads to an exact spectrum 
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FIG. 3: Second order structure functions S% ± (£) (solid) and FIG. 4: S£ x (solid) and (dot) compensated by £ 1/2 . The 
S2\\W ( dot )- Several slopes are given as a reference. The inset shows the same quant ities compensated by £ 2/3 . 
inset shows the ratio S%± (£) / S^u {£) in linear coordinates. 



~ fcj 2 /(fc||), where ku and k± are the amplitudes of the 
wave vectors in the directions parallel and perpendicular 
to a uniform magnetic field Bo- It has long been argued 
that at high Reynolds number, leading to a large scale 
separation, the large-scale magnetic field Bls plays the 
same dynamical role as Bo- As a result, we can expect 
the small scales to be anisotropic with respect to Bls- 
However, since Bls evolves in time with a characteristic 
correlation time tls and has a curvature proportional to 
the Taylor scale Am, in the absence of a uniform field 
and for isotropic initial conditions, one may assume that 
isotropy is still preserved in the large and intermediate 
scales (£ > Xm). 



To illustrate this, we follow the procedure given in [14 1. 
A local mean field Bls ( x ) is computed as the average of 
the magnetic field b in a box of size Lm with center on 
x, such as Bls = 3 )box(L. M xV Then, unit vectors 1|| and 
1_l (|| and _!_ with respect to Bls) are defined, and the 
two following longitudinal magnetic structure functions 



S 2 b iW = 5 2 b («i), 



(3) 



are computed for all points in the vicinity of x and for 
£ e fink'^, Lm/2]- The process is repeated for sev- 
eral values of x, and S 2± and are averaged over all 
the boxes (the results shown are the average over pa 10 7 
points). Figure [3] gives the resulting S%x and Shi, as 
well as their ratio. In the limit £ — > 0, this ratio can be 
associated with one of the so-called Shebalin angles [l^ 



tan 2 (0) = 2hm%^ 



J 2\\ 



(£)' 



(4) 



which measures the anisotropy of b with respect to the 
local mean field. 

We see that scales smaller than £ ps 0.4 are anisotropic, 
in agreement with \~\~\ . l?2_lW > S^(£) in this range. 



However, isotropy obtains in an intermediate range of 
scales corresponding roughly to the —3/2 range displayed 
in Fig. [TJ For very small increments {£ < 0.02 pa 2Trkp 1 ), 
both structure functions follow ~ £ 2 scaling, as expected 
for smooth fields in the dissipative range. At scales larger 
than the dissipation scale, but smaller than £ pa 0.4, a 
range with S^x ~ ^ is observed, consistent with a fcj 2 
scaling in the energy spectrum as predicted by WT the- 
ory [13[ ■ In the range £ ~ [0.4, 1.21 , where fluctuations are 



roughly isotropic, S^x ~ ~ Figure H] shows the 
structure functions compensated by £ x l 2 (for IK scaling) 
and £ 2 / 3 (for K41). The - l 1 ' 2 scaling in this range is 
again consistent with the isotropic energy spectrum and 
scaling of the isotropic (p = 3 and 4) structure functions 
discussed previously. 

The recovery of isotropy observed at intermediate 
scales is not obvious. Note that tan 2 (6') > 1, and for very 
small £, we have S 2 ± ~ ~ £ 2 . If Bls is strong, and 
there is enough scale separation for the local mean field to 
look as a uniform field for the small scales, WT theory is 
consistent with S 2 \\(£) ~ i 1/2 - Now assume S 2 x(£) ~ i a - 
For a > 1/2, the structure functions satisfy S2x/ S 2 \\ > 1 
at all scales. Thus, the recovery of isotropy at intermedi- 
ate scales should not be associated with the small scale 
fluctuations but rather with the fluctuations of the large 
scale magnetic field, and thus may not always take place. 
This is consistent with results from MHD with a dc mag- 
netic field: in this case, the fields are anisotropic at all 
scales [lij]. It may mean that in MHD the energy spec- 
trum is not universal and depends on some property of 
the large scale field Q. Extensive high resolution para- 
metric studies are needed to determine conclusively what 
is happening. 

Let Am be the largest scale for the fluctuations to see 
the mean magnetic field as a uniform field; for £ pa \mi 
the curvature of the large scale magnetic field cannot be 
neglected, and fluctuations become more isotropic. At 
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scales larger than Am, we thus have isotropic fluctua- 
tions with a local mean magnetic field Bls- Assuming 



equipartition and constant energy flux e 



(where tnl^ ~ £/vi is the nonlinear turnover time, and 
t~a,£ ~ £/Bls is the Alfven crossing time), we obtain 
bi ~ V£ ~ I 1 ! 2 and the IK energy spectrum E(k) ~ 
£;~ 3 / 2 . At scales smaller than Am, the fluctuations are 
anisotropic and several phenomenologies have been put 
forward 11, 13, Assuming vi and 6; are mostly 

in the direction perpendicular to Bls, Tnl.; ~ 
The anisotropic extension of IK takes the Alfven time as 



-1/2 



and the 



t~a,i ~ ^|| /-Bls, which results in ~ 6; - <._!_<., 
energy spectrum E(k) ~ fcj 2 fc^ 2 . In the context of this 
run, where the correlation length of the large scale mag- 
netic field is Lm, we can assume ta,z ~ Lm/Bls, which 
leads to the same ~ t±_ and ~ kj 2 scaling. Note that 
this _L spectrum is what WT theory predicts, and that 
for the isotropic case (fcj_ ~ fey), E(k) ~ kj_ 2 k^ 2 turns 
into the IK spectrum. In this light, it is not surprising 
that the IK spectrum at large scales be followed, in the 
same dynamical vein, by its wave turbulence anisotropic 
counterpart at small scales. 

Another spectrum for MHD turbulence has been advo- 
cated in [l6| , whereby the anisotropy of the flow induces a 
Kolmogorov spectrum. In the anisotropic range, this im- 
plies (Sgj, ~ t 2 ^ ■ Another anisotropic phenomenological 
model based on field alignment [ll| implies S\±_ ~ I 1 ! 2 . 
Such scalings are not observed in our computation, and 
the results in Fig. [4] suggest S\±_ ~ I is a better fit in the 
anisotropic range. However, other scaling laws cannot be 
ruled out completely for the following reason. It is well 
known that the assumption leading to WT (namely, that 
the characteristic time of the wave be shorter than the 
time for nonlinear exchange of energy between eddies) 
is non-uniform in scale, and breaks down at small scales 
when the turn-over time has become sufficiently small. It 
is possible that we have not achieved yet sufficient scale 
separation for this to occur, and simulations at higher 
Reynolds number are needed in order to reach a better 
and more detailed understanding of such flows. 

Since we also know that MHD turbulence displays non- 
local transfer of energy [l7j], leading to a direct connec- 
tion between the large and the small scales, a break-down 
of universality can be a candidate to explain the variety 
of solutions reported in the literature. The parameters 
that can influence such a break-down are many and will 
require further study. In this context, it is useful to com- 
pare the present results with simulations of forced MHD 
turbulence. In our simulations we can associate the scale 
where the transition from the isotropic to anisotropic 
field takes place with Am, and the correlation time of 
the local mean field tls with Lm/B^q. In a forced sim- 
ulation, the correlation time of the large scale magnetic 
field can be associated with the correlation of the exter- 



nal forcing rp. If rp <C ta,£, then the local mean field 
changes faster than the Alfven wave crossing time at the 
scale £, and we should not expect anisotropy to develop 
at that scale, or wave interactions to be relevant. The 
scale where the transition takes place would change ac- 
cordingly. Such a dependence would be consistent with 
results obtained in forced RMHD simulations 0] and re- 
cent MHD simulations [Tj |. 
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